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ABSTRACT 


Through  the  use  of  asymptotic  and  perturbation  methods, 
this  thesis  presents  a  theoretical  study  of  the  flow  of  a 
stratified  fluid  over  variable  topography  as  a  model  of  the 
resonant  interaction  of  a  s\ibmarine's  wake  with  a  stratified 
fluid.  Such  resonant  interactions  may  be  able  to  produce 
significant  upstreaun  disturbances.  The  long  time  solution 
obtained  in  our  model  exhibits  growth  in  time  for  the  resonant 
case,  indicating  that  perhaps  nonlinear  effects,  balanced  by 
dispersion,  could  cause  significant  upstreaun  disturbances. 
These  disturbances  could  allow  the  possibility  of  detecting  a 
submarine  by  the  internal  waves  its  wake  generates  in  the 
resonant  case . 
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I .  ISTEQDUCTIC^ 


The  nonacoustic  detection  of  submarines  is  an  area  of 
intensive  and  well  funded  research.  Although,  it  is  not  clear 
whether  devices  based  on  nonacoustic  means  of  detection  will 
become  as  useful  as  those  based  on  acoustic  means,  the 
military  in^ortance  of  such  devices  demands  pursuing  such 
avenues  of  technological  research.  One  possible  nonacoustic 
means  of  detection  is  to  attempt  to  detect  the  submarine 
optically  through  the  use  of  a  laser  by  examining  the  way  the 
submarine  reflects  or  absorbs  blue-green  light.  Another  more 
promising  method  is  to  try  to  detect  the  submarine  by  the  heat 
it  emits.  However,  an  even  more  promising  technique  is  the 
detection  of  submarines  by  the  internal  waves  they  generate. 

In  the  ocean  a  stable  density  stratification  generally 
exists.  As  the  submarine  moves  through  the  density- stratified 
ocean,  it  leaves  behind  a  wake  capable  of  generating  internal 
waves.  These  internal  waves  result  from  the  stedsilizing 
effect  of  the  stratification:  as  a  water  particle  is 
displaced,  it  experiences  a  buoyant  restoring  force  in  the 
opposite  direction  from  which  it  was  displaced. 

Many  theoretical  and  experimental  studies  have  been 
conducted  on  the  subject  of  the  flow  of  a  stratified  fluid 
over  a  localized  topography.  For  the  two-dimensional  flow  of 
an  inviscid,  incompressible  stratified  fluid,  most  studies 


have  been,  confined  to  steady  flow  or  linearized  theory. 
Earlier  studies  concerned  themselves  only  with  a  description 
of  the  downstream  stationary  lee-wave  field.  A  question  that 
has  generated  some  controversy  is  the  phenomenon  of  upstream 
influence.  To  what  extent  do  disturbances  generated  by  the 
topography  through  transient  and  nonlinear  processes  propagate 
upstream? 

McIntyre  (1972)  showed  that  the  only  significant  upstream 
disturbances  were  weak  second  order  long-wave  motions 
generated  by  nonlinear  interactions  in  the  lee-wave.  McIntyre 
(1972)  extended  the  earlier  theories  of  Benjamin  (1970)  and 
Ready  (1971)  by  performing  an  expansion  using  as  the  small 
parameter  the  ratio  of  the  height  of  the  topography  to  the 
total  fluid  depth.  Baines  (1977,1979),  on  the  other  hand,  in 
experimentally  excunining  the  flow  of  a  continuously  stratified 
fluid  over  topography,  found  upstream  disturbances  that  were 
first  order  in  the  topography  height.  These  upstream 
disturbances  were  not  only  generated  by  nonlinear  processes 
over  the  topography  but  also  travelled  at  the  long- wave  speed. 
Baines  also  observed  that  these  upstream  disturbances  were 
very  strong  when  the  fluid  approached  resonance.  Here, 
resonance  is  defined  as  the  near  coincidence  of  the  basic  flow 
speed  and  one  of  the  free  long-wave  speeds.  Grimshaw  and 
Smyth  (1986)  presented  a  theoretical  study  of  stratified  fluid 
flow  over  localized  topography  for  the  two-dimensional  case 
when  the  basic  flow  is  near  resonance.  One  interesting 
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feature  of  their  results  is  the  appearance  of  significant 
upstream  disturbances  as  observed  by  Baines  {1977,1979}. 
Grimshaw  and  Sayth  showed  that  solitary  waves  were  formed  when 
the  localized  forcing  over  the  topography  produced  a  positive 
disturbance  upstream,  and  an  oscillatory  wavetrain  was  formed 
when  the  disturbance  had  negative  polarity. 

Motivated  by  these  results,  this  thesis  presents  a 
theoretical  study  of  the  flow  of  a  stratified  fluid  over- 
localized  topography  for  the  three  dimensional  case. 
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ZI.  PTOBm  F0SllDZAZ>Z0i9r 


K.  DIMBISlCaaL  PKOBLSK 

We  consider  the  three-dimensional  flow  of  an  inviscid, 
inccHi^ressible  stratified  fluid.  Figure  1  (see  page  5) 
provides  a  sketch  of  the  coordinate  system.  He  assun^  that 
the  basic  state  has  a  constant  horizontal  velocity,  V,  frcan 
left  to  right  with  equilibri\ain  pressure  and  density 
respectively  as  p^tz)  and  pgiz)  .  Here  and  throughout  this 
paper  barred  variables  represent  dimensional  variables.  We 
introduce  f  as  the  vertical  particle  displacement  from  some 
initial  position  at  time  t  =  0  to  the  position  at  time  t  •  T. 
Using  the  Lagrangian  represents  ion  f(x,y,z,T)  -  Zf  -  z„,  the 
density,  p(x,y,z,t),  beccanes  po(z-f  (x,y,  z,  t) )  .  Therefore,  the 
equations  of  motion  are 


Uj  +  V-  +  =  0, 

P„(z-f)  —  +  ft  =  0, 

m 

^ft  =  0, 


^  p^ 


z-f)g  +  p  =  0, 


(2.1) 

(2.2) 

(2.3) 

(2.4) 
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where 


(2.5) 


Dt 


±  . 

art 


a  a  *”  d 

( V+u) -  ♦  V -  +  w - 

Bx  ay  dz 


(2.6) 


Here,  u,  v,  and  w  are  the  velocity  cotr^yonents  relative  to  the 
basic  state  which  has  velocity  (V,0,0)  with  respect  to  a 
right-handed  Cartesian  coordinate  system. 

The  bottom  topography  is  given  by  z  -  f(x,y,t),  where 
initially  f(x,y,0)  ••  0.  For  large  times  the  topography 

attains  a  steady  value  g(x,y}  and  hence  we  assume  that 

f(x,y,  t)  -  ir(x,y)  as  t ®.  (2.7) 


•  • 


We  also  assume  that  at  all  times  the  topography  is  localized 
where  f  {x,y,t)  -•  0  as  }?+y  -»  ».  The  kinematic  bottom  boundary 
condition  is  then 


f  =  f  at  z  =  f. 


(2.8) 
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*  •  •  •  •  •  •  . m...  •  • 


We  are  assuming  that  this  "topography*  is  created  by  the 
submarine's  %»ake,  which  when  viewed  frtsn  the  moving  sul:»narine, 
actually  has  a  steady  shape.  In  this  thesis,  we  investigate 
the  interaction  of  the  wake  with  internal  waves.  KSllen 
(1987)  has  investigated  this  type  of  interaction  in  a  linear, 
steady  state  nuDdel  and  comments  in  detail  on  the  validity  of 
assuming  that  the  wake  acts  as  rigid  topography.  Furthermore, 
from  experimental  results,  Lin  and  Pao  (1979)  show  the  mean 
flow  occurs  at  the  wake  boundaries  implying  that  part  of  the 
flow  travels  around  rather  than  through  the  wake.  It  is  this 
flow  around  the  wake  that  generates  the  internal  waves.  Our 
model,  in  fact,  does  not  depend  critically  on  the  assumed  form 
of  "topography"  the  submarine  creates  and  our  purpose  here  is 
to  investigate  a  time  dependent  model. 

The  upper  boundary,  assumed  to  be  a  free  surface,  is 
located  at  z  «  H  +  J)(x,y,t),  where  H  is  the  vertical  distance 
from  the  surface  to  the  submarine  and  ^(x,y,lD)  is  the  free 
surface  displacement.  One  must  remember  to  include  not  only 
the  kinematic  boundary  condition  at  the  upper  boundary  but  the 
dynamic  boundary  condition  as  well.  The  upper  boundary 
conditions  take  the  form: 


•  • 


f  =  »j  at  z  =  H  *  i\. 


(2.9a) 


Po(H+ij)  +  p,(H+»?)  =  p„(H)  at  z  =  H  + 


(2.9b) 
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Here,  p,  is  the  perturbation  to  the  pressure  p*  of  the  basic 
state. 

B.  DIMBKSIOiaU.  SCAX.IHG 

The  solution  of  practical  problems  in  fluid  mechanics 
requires  both  theoretical  developments  and  estperimental 
results.  By  grouping  significant  dimensional  quantities  into 
dimensionless  parameters  it  is  possible  to  present  results  in 
a  compact  form  which  is  applicable  to  all  >3imilar  situations. 

In  the  problem,  the  characteristic  horizontal  length  scale 
L  is  determined  by  the  ratio  U/N,  where  U  is  a  characteristic 
horizontal  velocity  and  N  is  the  Brunt -VS,is&lS  frequency  at 
the  location  of  the  submarine.  We  shall  take  the 
characteristic  horizontal  velocity  U  to  be  given  by  N,K,  (U  » 
N,H)  ,  where  N,  is  the  maximum  Brunt-VSisSlS  frequency  and  H, 
the  characteristic  vertical  length,  is  the  depth  of  the 
submarine.  Typically  N,  is  10’^  sec"'  and  H  is  200  meters  (m)  , 
so  U  is  2  m/sec.  This  is  a  typical  order  of  magnitude 
patrolling  speed  for  a  submarine  at  that  depth.  Of  course, 
0(1)  variations  in  N,  and  H  can  be  made  without  changing  our 
analysis. 

From  the  mass  conservation  equation  and  with  the  inqposed 
horizontal  length  and  velocity  scales,  we  see  that  the 
vertical  velocity  scale,  W,  is  UH/L.  The  horizontal  length 
and  velocity  scales  are  also  used  in  determining  the 
characteristic  time  scale,  T,  where  T  is  defined  by  T  -  L/U  » 

B 


l/N.  Lastly,  the  characteristic  pressure  scale  is  defined  as 
P  =  Pi3H  where  is  a  characteristic  density.  These  scales 
define  the  parauneter  /S  =  N,^/g,  which  is  a  dimensionless 
measure  of  the  ratio  of  the  buoyancy  force  to  the 
gravitational  force. 

C.  HCaiDIKBMrSICmLIZATiaN 

By  using  the  various  dimensional  scales  introduced  in 
section  B,  we  shall  nondimensionalize  the  equations  of  motion 
(2. 1-2. 6)  along  with  the  boundary  conditions  (2. 8-2. 9b). 
However,  before  deriving  the  nondimensional  equations  of 
motion  and  boundary  conditions,  we  will  consider  several 
related  dimensional  equations  that  will  assist  in  simplifying 
the  nondimensional  derivation  of  equations  (2. 1-2. 6)  and 
boundary  conditions  (2. 8-2. 9b).  We  first  introduce  the 
following  nondimensional  variables: 


% 


•  i 


N  = 


(2.10) 


To  nondimensionalize  the  hydrostatic  equation,  we  begin  with 


cy=>« 

dS 


=  -9Po- 


Using  (2.10)  yields 


(2.11) 


# 


Po(2)  =  -/>»• 


(2.12) 


-  • 


The  dimensional  Brunt -VSisSia  frequency  is  defined  as 


# 


(2.13) 


where  dp^/dz  <  0.  Using  (2.10),  we  obtain 


1  dPo 


(2.14) 


but  |5  -  N,*H/g,  therefore  (2.14)  sinplifies  to 


=  - 


~3F’ 


(2.15) 
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In  nondimensionalizing  the  Lagrangian  representation  of  the 
particle's  vertical  displacement,  we  start  with 

«„  =  z  -  f  (x,y, "i,  t)  .  (2.16) 

Using  (2.10)  yields 

z„  =  z  -  f  (x,y,  2,  t)  .  (2.17) 

The  last  equation  needed  before  considering  the  equations  of 
motion  and  boundary  conditions  is  the  dimensional  pressure 
equation: 

p(x,y,z,t)  =  p,gHp„{2)  *  fiygll&q(x,y,z,t)  (2.18) 

where  pp  »  p,gH)3q (x,y,  z,  t)  and  &  is  defined  after  (2.14). 
Thus,  q(x,y,2,t)  is  the  dimensionless  perturbation  to  the 
hydrostatic  pressure  p,,. 

Before  deriving  the  nondimens ional  equations  of  motion  and 
their  associated  boundary  conditions,  we  first  introduce  the 
following  nondimens ional  variables: 


X  = 

X 

L  ' 

y  =  i 

^  L  ' 

V  = 

V 

U  ' 

-1  ■ 

“=1 ' 

Pi 

t  = 

t 

T  ' 

r  =  i  , 

■ 

(2.19) 


"  =  1 
H  ^  ' 


where  L,  H,  U,  W,  p,,  and  T  are  the  characteristic  dimensional 
scales  defined  in  section  B  and  h  is  a  characteristic 
anqplitude  of  Che  topography. 

We  begin  with  the  mass  conservation  equation  (2.1)  .  Using 
(2.19)  and  remembering  that  W  *  UH/L,  (2.1)  becomes 


u,  +  V,  +  Wj  =  0. 


(2.20) 


To  assist  in  deriving  the  nondimens ional  form  of  the  x- 
coii5)onent  of  the  momentum  equations,  given  by  (2.2),  we  shall 
now  derive  the  nondimens ional  form  for  the  convective 
derivative  D/Dt: 


D  d  iTr  d  —  3  —  3 

—  =  +  (V+u)—  +  V—  +  w— , 

Dt  3'C  3x  3y  dz 


(2.21) 


Using  (2.19),  (2.21)  becomes 
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but  T 


£>  ^  1  a 
m.  T  Ti 


+ 


M  ( v+u) 

jj 


d 

'7k 


+ 


-  L/U,  therefore  (2.22) 


D 

IRE 


0  d  0  d 
L  ay  L  0z 


simplifies  to 


1  D 
T  Dt' 


where 


(2.22) 


(2.23) 


Dt  Tt 


(V+u) 


d 

■33i 


+  V 


d 

■5y 


+  w. 


Sk- 


(2.24) 


Substituting  (2.17),  (2.18),  and  (2.23)  into  (2.2)  and  using 

(2.19),  we  obtain 


9 


9 


m 


# 


PlP„(2-f) 


0  0U 
T  Dt 


+ 


p,gH0 

L 


<gt  =  0  or. 


(2.25) 


Po 


+  q,  =  0. 


A  similar  analysis  for  the  y- component  of  the  momentum 
equations  (2.3)  yields 

« 
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+  g,  =  0. 


(2.26) 


To  derive  the  nondimensional  equation  for  the  z-ccm^onent  of 
the  moBientum  equations  (2.4),  we  begin  with 


Pp(z-f)  ~  +  +  ft  =  0. 


(2.27) 


Substituting  (2.17),  (2.18),  and  (2.23)  into  (2.27)  and  using 
(2.19),  we  obtain 


+  P„(2-f)Sr  +  9P«  +  90<lz  =  0- 


(2.28) 


Using  T  =  L/N,H  and  W  -=  N,hVl,  (2.28)  simplifies  to 


P<,{2-i')3'+  9P«  +  SflSg  =  0.  (2.29) 


Using  (2.12)  and  the  definition  ot  0  =  N,*H/g,  (2.29) 

sin^lifies  to 


a 


Jj 


(2*35} 


For  the  free  surface  kinematic  boundary  condition. 


# 


# 


^  =  1}  at  z  =  H+i}. 


(2.36) 


'  • 


We  use  (2.19)  to  obtain 


f  =  at  z  =  l  +  ^.  (2.37) 

Next,  we  nondimensionalize  the  free  surface  dynamic  boundary 
condition  in  (2.9b); 


pJH*ri)  +  Pp(H+i7)  =  p„(H)  at  z  =  H  +  II.  (2.38) 

Using  (2.19)  and  (2.18),  (2.38)  becomes 

p„(l-t-if)  +  0q  =  p„(l)  at  z  =  1  +  If.  (2.39)  ^ 

To  recapitulate  our  nondimens ional  versions  of  equations 
(2. 1-2. 6)  amd  the  surface  boundary  conditions  (2. 8-2. 9b),  the  • 

nondimens ional  equations  of  motion  are 
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V,  Wj  =  0, 

(2.40) 

*3,  =  0. 

(2.41) 

• 

P.u-n^  ^g,  =  o. 

(2.42) 

-*■  ■J(Po(*-r)-p»(«)}  +  =  0. 

(2.43) 

• 

w= 

ct' 

(2.44) 

• 

where  e  =■  H/L  and 

• 

D  3  /  rr  »  3  3  3 

-  =  +  (V+u)  +  v-J_  + 

Dt  7t  lx  ly  Iz' 

(2.45) 

• 

with  boundary  conditions 

i 

f  =  af  at  z  =  orf. 

(2.46) 

• 

f  =  il  at  z  =  l+i>. 

(2.47a) 

p„{l*V)  +  =  P,(l)  at  z  =  1  +  ty. 

(2.47b) 

• 

where  a  «  h/H. 

• 

An  exact  analytical  solution  of  these 

nonlinear  partial 

differential  equations  is  not  available. 

However,  we  shall 
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• 

obtain  approximate  solutions  for  the  case  when  e  defined  by 
(2.31)  and  a  defined  by  (2.35)  are  both  small  positive 
constants.  Since  e  «  H/L  »  N/N,,  the  condition  e  «  1  is 
satisfied  when  the  density  changes  rapidly  across  a  thin 
thermocline.  However,  there  is  scane  ambiguity  in  the 
definition  of  the  thermocline.  He  will  consider  the 
thermocline  to  be  the  region  where  density  variations  are 
determined  chiefly  by  ten^jeratxire  variations.  This  means  that 
the  siiJxsarine  is  usually  sulanerged  to  a  depth  below  the 
thermocline.  Is  this  a  valid  assumption?  To  avoid  acoustic 
detection  the  sulxnarine  submerges  to  depths  below  the  mixed 
layer  or  layer  depth.  Depending  on  the  thickness  of  the 
thermocline  the  sulanarine  cem  be  either  in  the  thermocline  or 
belcw  it.  We  shall  consider  only  the  case  for  which  the 
submarine  is  below  the  thermocline  in  which  the  pareuneter  a, 
the  ratio  of  the  characteristic  anplitude  of  the  topography  to 
the  siiixnarine ' s  depth  defined  by  (2.35)  ,  is  also  small  in  this 


situation. 


111.  TBB  JtSTKPTOTlC  B3CPAXfSI(»I  OF  TBS  DBFEODBHT  VS&XABLBS 


To  study  the  behavior  of  the  internal  waves  excited  by  a 
small  amplitude  disturbance  in  the  presence  of  a  uniform 
oncoming  flow,  as  when  a  suixnarine  moves  through  the  water  at 
constant  velocity  V,  we  perform  an  asynptotic  expansion  on  the 
variables  f,  q,  u,  v,  and  ij  in  terms  of  the  small  parameter  <*- 
This  is  reasonable  since  the  small  parameter  a  enters  the 
bottom  bo\indary  condition  in  an  analytic  way  and  since  the 
solutions  arise  by  analytic  processes,  it  is  natural  to  expect 
the  solution  to  exhibit  an  analytic  dependence  on  the 
parameter  a.  Since  the  disturbance  is  0(a),  we  shall 
initially  suppose  that  the  response  is  also  0(a),  Therefore, 
we  introduce  the  following  asynptotic  expansions 


f  =  af^(x,y,  z,  t)  +  +  . . . 

q  =  aq^ix.y,  z,  t)  +  opq,  +  . . . 

u  =  au„(x,y,2,t)  +  Q(^u,  +  ...  (3.1) 

V  =  av^(x,y,  2,t)  +  a^v,  +  — 

ij  =  (x,  y.  2,  t)  +  +  . . . 

in  which  we  consider  the  limit  process  as  a  0. 

Substituting  (3.1)  into  the  continuity  equation  (2.40),  we 


obtain 


•v5i» 


-^[auj  ♦  +  0{a»)  +  ...  =  0.  (3.2) 

Dividing  by  a  and  taking  the  limit  as  a  -*  0,  (3.2)  simplifies 
to 

+  w;,  =  0,  (3.3) 

where 


and 


w. 


(3.4) 


-fft  Sx' 


(3.5) 


To  show  that  D^Dot  is  the  0(1)  part  of  the  convective 
derivative  given  in  (2.45),  we  collect  only  the  0(1)  terms  in 
the  expansion  for  D/Dt  given  in  (3.6)  below: 


-5-  =  «  . 

Dt  7t 


(V  +  au„) 


d 


+  QfV, 


d 

■3y 


O(a^)  +  . . . 


(3.6) 
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Eubstituting  (3.1)  into  the  x-cc»nponent  of  the  momentuin 
equations  (2,41),  we  obtain 

(ofuj  +  aq„  +  O(cK^)  +  . .  .  =  0.  (3.7) 

Expanding  p„(Z“f)  in  a  Taylor  Series  and  using  (3.1)  yields 

P«(z-f)  =  P<,(z)  "  Po(z)  [afj  +  0(Q!^)  + -  (3.8) 

Substituting  (3.8)  into  (3.7),  we  obtain 

Otp^iz)  +  aiq^  +  0(0?^)  +  .  .  .  =  0.  (3.9) 

Dividing  by  a  and  ta)ting  the  limit  as  a  -*  0,  (3.9)  becomes 

=  C.  (3.10) 

A  similar  analysis  for  the  y- component  of  the  momentum 
equations  (2.42)  yields 
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(3.11) 


*a,  =  0. 


For  the  z-ccffiponeiJt  of  the  momenttim  equations  (2.<3),  we 
substitute  (3.1)  into  (2.43)  to  obtain 


^  [p  (z-f)  -p^(z)]  +  aq^  +  0(Qf*)  +  .. 

n  t®  P 


=  0.  (3.12) 


Expanding  p„(z-f)  in  a  Taylor  Series  as  in  (3.8)  and 
substituting  into  (3.12),  we  obtain 


t* 

-  ■pPa<2)  +  a<l„  +  0(«*)  +  •  •  •  =  O' 


(3.13) 


•  • 


Dividing  by  a  and  talcing  the  limit  as  a  -*  0,  (3.13)  simplifies 
to 


fo  +  =  0. 
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(3.14) 


Notice  that  the  e®  term  has  been  neglected.  As  previously 
discussed,  in  order  for  our  theory  to  worlc  we  require  «  «  i. 
Therefore,  we  see  the  e®  term  is  small.  The  fact  that  the  £® 
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41- 


term  is  small  by  itself  is  not  justification  in  dropping  this 
term  frcan  (3.13).  In  order  to  justify  the  dropping  of  the  £“ 
term,  we  shall  assume  that  =  o(a)  as  a  -*  0,  or  that  <* 
0(1)  as  a  -*  0.  Now  as  one  takes  the  limit  a  -*  0  in  (3.13)  ,  we 
clearly  see  that  0.  In  continuing  our  derivation  of  the 
z - cort^ionent  of  the  momentum  equations,  we  sxibstitute  (2.15)  in 
for  dp<,/dz,  so  that  (3.14)  becomes 

(3.15) 

Next,  we  consider  the  surface  boundary  conditions.  We 
begin  by  substituting  (3.1)  into  the  kinematic  bottom  boundary 
condition  (2.46)  to  obtain 


-  af  +  0{o^)  +  . . .  at  z  =  af.  (3.16) 

Taking  the  limit  as  a  -*  0,  we  find 

=  f  at  z  =  0.  (3.17) 

Note  that  the  0(a)  bottom  boundary  condition  is  applied  at  z 
»  0.  Siibstituting  (3.1)  into  the  free  surface  kinematic 
boundary  condition  {2.47a)  results  in 


affl  =  <x%  +  0(0^)  +  . . 


at 


*  =  1  +  aij„. 


{3.18) 


:  • 


Taking  the  limit  as  a  -*  0,  (3.18)  becomes 


“  « 


i;  =  n,  at  2  =  1.  (3.19) 

;;  • 

Again  the  0(a)  free  surface  condition  is  applied  at  z  «  i. 

The  last  boundary  condition  to  consider  is  the  free  ^ 

surface  dynamic  boundary  condition  (2.47b).  Siibstituting 
(3.1)  into  (2.47b),  we  obtain 


p^d+ij)  +  +  0(0^)  +  .  .  .  =  p„(l)  .  (3.20) 

Expanding  p^d+ij)  in  a  Taylor  Series  and  using  (3.1)  gives 

p„d+tj)  =  p„d)  +  ap!,(z)ii„  *  O(o^)  +  ...  (3.21) 

Substituting  (3.21)  into  (3.20),  we  find 

p^d)  +  ap' (l)ii„  +  ajSqr,  +  O(a^)  +  ...  =p„(i).  (3.22) 
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Canceling  the  Po(l)  terms,  dividing  by  a,  and  taking  the  limit 
as  Of  0,  (3.22)  beccmies 


Pa  +  /Sq;,  =  0  at  s  =  l  (3.23) 

but  using  (2.12)  sin5>lifies  the  free  surface  dynamic  boundary 
condition  to 


PQo  =  PcVo  at  z  =  1.  (3.24) 

In  summary,  we  see  from  sxibstituting  the  asyn^totic 
expansion  (3.1)  into  (2.40-2.45),  (2.46),  and  {2.47a,b),  that 
the  0(a)  problem  is; 


•  i 


=  0, 

(3.25) 

D.u. 

-ffr 

(3.26) 

D.v 

^  0, 

(3.27) 

<I«  +  Pj^^o  =  0, 

(3.28) 
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(3.29) 


_ 

_ 


where 


5  +  V  « 
ID^  7t 


(3.30) 


with  boundary  conditions 


to  -  f 

at 

o 

II 

(3.31) 

to  =  Vo 

at 

z  =  1, 

(3.32a} 

=  PoVo 

at 

Z  =  1. 

(3.32b) 

Combining  (3.32a)  and  (3.32b}  gives  the  single  boundary 
condition  at  the  free  surface: 


P<1,  =  Po^o  at  z  =  1. 


(3.33) 
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nr.  DKvxLOPMBsrr  of  the  Two-DxicBaiszoauiL  wavs  bquxtxow 

A.  BI6ENVALUB  FROBLSM 

In  order  to  solve  (3.25-28),  (3.31),  and  (3.33),  we  look 
for  solutions  of  the  form 


Ug  =  Ae‘ « -  or )  (z) 

«(z) 

g,  =  p„  *Mz) 


(4.1) 


where  k  »  (K,L)  is  the  horizontal  wave  vector,  0  is  the 
frequency,  and  c  is  the  long-wave  phase  speed  relative  to  the 
basic  state  at  rest.  Substituting  (4.1)  into  (3.25),  (3.26), 
and  (3.27)  and  solving  the  system  of  linear  equations,  we  find 
the-  dispersion  relationship  to  be 


Q  =  cjK(  where  |itj^0.  (4.2) 

Furthermore,  substituting  (4.1)  into  (3.28),  (3.31),  and 

(3.33)  results  in 
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(4.3a] 


=  0, 


with  boundary  conditions 


4t(z)  =0 

at 

z  =  0, 

(4.3b) 

0(z)  =  (z) 

at 

z  =  1. 

{4.3c) 

The  equations  (4.3a-c)  formally  pose  am  eigenvalue  problem 
where  c'^  is  the  eigenvalue  and  ^  the  eigenfunction.  Note  that 
the  dispersion  relationship  in  (4.2)  was  obtained  for  wave 
motion  in  a  stratified  fluid  at  rest  where  V  «  0  and  D„/Ti„t  =« 
d/dt.  In  the  problem  we  are  considering,  the  dispersion 
relationship  would  differ  from  that  found  in  (4.2)  but  the 
eigenvalue  problem  would  remain  the  saune.  Hence,  our  purpose 
here  is  to  show  the  existence  emd  character  of  the  solutions 
that  define  the  eigenvalue  problem.  The  solutions  have  the 
horizontal  character  of  plane  waves  travelling  in  the  (x,y) 
plane  with  a  vertical  structure  defined  by  the  eigenfunctions, 
^(z)  .  We  shall  see  that  the  general  solution  to  (3.25-3.28), 
(3.31),  (3.33)  can  be  obtained  in  terms  of  a  superposition  of 
normal  modes  whose  depth  dependence  is  determined  by  the 
0, (z) ,  the  eigenfunctions  of  the  following  eigenvalue  problem: 


28 


O 

II 

+ 

(4.4a} 

C, 

o 

II 

14 

at  z  =  0, 

(4.4jb) 

#,(z)  =  0cf4>„{z) 

at  z  =  1. 

(4.4c) 

Here  c,'^  is  the  eigenvalue,  where  c,  is  the  long-wave  phase 
speed  relative  to  the  basic  state  at  rest. 

B.  OBTBOGCMALITY 

The  eigenfunctions  4>,{z)  satisfy  the  orthogonality 
condition 


A  .  (4.5) 

where  I,  is  defined  by  the  left  side  of  (4.5)  for  r  »  s.  To 
prove  (4.5) ,  we  let  X,  and  X,  be  eigenvalues  with  corresponding 
eigenfunctions  and  4>r  where  X,  -  c,'^  and  X,  -  c,'^.  The 
differential  equations  satisfied  by  the  eigenfunctions  are 


)'  +  =  o» 

(4.6a) 

(pA)'  \Pot^*r  =  0. 

(4.6b) 
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In  addition,  both  4^,  and  4>t  satisfy  the  same  set  of  homogeneous 
boundary  conditions  defined  in  <4.4b,c).  S&iltiplying  (4.6a) 
by  and  (4,6b)  by  yields 


K(Po<t>^  ) '4>r  +  KKPa^^<t>. =  0 ,  (4 . 7a) 

\  p,ifi4>,4>r  =  0.  (4 . 7b) 

Subtracting  (4.7b)  from  (4.7a)  results  in 

\  )'4>r-K  Mn  )'4>,  =  0.  (4.8) 

Through  application  of  the  product  rule 

K^PcK)<t>r  -\UP.K<t>ry  -Pc*n<f>n^'  (4.9) 

(4.8)  becomes 

K’^Po<l>K<t>ry  -  KPoK’t>rz  -  K^p0<t>n<l^>y  *  \P0KK  =  0-  (4.10) 


Integrating  (4.10)  yields 


(4,11) 


Evaluating  the  right  hand  side  of  (4.11)  and  remembering  ^,(1) 
-  .  ^,(1)  -  /3V‘«„{1).  and  ^,(0)  -  ^.(0)  -  0,  (4.11) 

beccm^s 


(4.12) 


If  X.  it  follows  that 


(4.13) 


In  other  words,  the  z  derivatives  of  the  eigenfunctions  (^i, 
and  0.)  corresponding  to  different  eigenvalues  (X,  and  X,)  are 
orthogonal  with  weight  p„.  Therefore,  the  normal  modes  do 
satisfy  the  orthogonality  condition  in  (4.5) . 

C.  TWO-DIMENSIONAL  WAVE  BQaATION 

Since  the  eigenfunctions,  ^,(z)  ,  form  a  complete  set,  we  • 

can  represent  u,,,  v,,  and  by  a  generalized  Fourier  series 
expansion.  From  (4.1)  we  )cnow  the  general  form  of  the 
solutions  and  hence  we  put  • 
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<4. 14a) 


Co  =  52A,(x,y.  t)^,(2) , 

0 

m 

<z) 5^  c|b, (x. y,  t> (z) ,  (4 . 14jb) 


Uo  =  D^ix.y.  t)  ^„(z) .  U .  14c) 

0 


Vo  =  -52£',(-?f»y.  t)^„<2) . 
0 


(4.14d) 


To  derive  a  relationship  between  A,  and  B,,  we  begin  with 


(4.15) 


We  take  the  derivative  of  (4.15)  with  respect  to  z,  multiply 

by  p„c,*0tt{z)  and  use  the  orthogonality  condition  in  (4.5)  to  * 

obtain 

• 

=  c,"/'poCo,«„dz-  (4.16) 

y  0 

Integrating  (4.16)  by  parts  then  gives 
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Using  (4.4a),  (4.17)  becomes 

X*A„  =  c-|po^«f  +  fyo(o^^*s 
0  *^0 

Now,  we  multiply  (4.14b)  by  and  f 
condition  (4.5),  we  obtain 

Integrating  (4.19)  by  parts,  we  obtain 

<3  Jo 

Using  (3.28),  (4.20)  becomes 

Combining  (4.18)  and  (4.21)  yields 
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orthogonality 


^IPo^o^srl  *  {4.22) 

0  0 

Substituting  (3.33)  and  {4.4c)  for  q,(l)  and  ^.(1)  respectively 
and  remembering  that  fo(0)  «  f(x,y,t),  (4.22)  sinqplifies  to 


m 


# 


B,  =  A,  +  F,. 


(4.23) 


where 


lgFg=  cjp^(0)^,j(0) /{x,y,  t)  .  (4.24) 

Note  that  P, (x,y,t)  is  proportional  to  the  bottom  topography 
f {x,y, t)  . 

Next,  we  substitute  (4.14a,c,d)  into  the  continuity 
equation  (3.25)  to  obtain 


A  Se 


0 


Dot 


-  0, 


(4.25) 


Multiplying  (4.25)  by  c.Vo^n  applying  the  orthogonality 

condition  (4.5)  then  gives 


•  • 
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az?,  as, 

I>^,t  ^x  dy 


(4.26) 


where  D(,/D,t  is  given  in  (3.30). 

To  derive  the  two-dimensional  wave  equation,  we  use  the  x 
and  y  coicqponents  of  the  mcsnentum  equations  (3.26)  and  (3.27). 
Substituting  (4.l4b,c)  and  (4.i4b,d)  into  (3.26)  and  (3.27) 
respectively,  we  obtain 


‘’O.t 

0 

0» 

0 

0 

^  cf  B,  <()„ 

0 

(4.27a) 


(4.27b) 


Multiplying  (4.27a)  and  (4.27b)  by  and  applying  the 

orthogonality  condition  (4.5)  then  yields 


+ 


.2  dB, 


=  0, 


(4.28a) 


jdB 

+  c’  *  =  0.  (4.28b) 

D„t  ^ 
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We  next  take  the  x  derivative  of  (4.28a)  and  the  y  derivative 
of  (4.28b)  and  add  the  two  results  to  obtain 

Substituting  (4.23)  and  (4.26)  into  (4.29)  finally  yields 


+  cJV^B.  =■■■■  0. 


(4.29) 


dIa, 


c/V^[A,  ^  F,]. 


(4.30) 


where 

(4.31) 

(4.32) 

Equation  (4.30)  is  the  inhomogeneous  two-dimensional  wave 
equation  whose  solution  will  be  obtained  in  Chapter  V. 

To  summarize,  we  see  that  substituting  (4.14a-d)  into 
(3.25-3.27)  implies 


•  « 
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Sf  —  Aj  +  Fj, 

^  ^  as, 

D„t  ~5x  ~Sy ' 


D„t'^ 


=  C,*r*[  A,  +  F,  ]  , 


(4.33a) 

(4.3313) 

(4.33c) 


where 


I,F,  =  c]pJ0)4>„mf(x,y,t),  (4.33d) 


and 


_  a 

D,t  Tt 


(4.33e) 


# 


3y2 


(4.33f) 


Here,  since  f  vanishes  at  t  -  0,  F,(x,y,0)  ■=  0  and  from  (2.7) 
we  have 


F,(x,y,  t)  "  G,{x,y) 


as  t  -*  oo, 


(4.34a) 


where 
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J,G,  =  c,V,(0)^„(0)g{x,y)  . 


(4.341}) 


;• 

The  initial  conditions  for  {4.l4a-d)  are  that  A,,  B,,  D,,  and 
E,  all  vanish  at  t  >  0. 
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V.  SOLUTZ(»r  TO  THB  TWO-DZMBlfSIOSnii  WAVS  BQUATI(»r 


A.  METHOD  OF  DBSCBMT 

To  solve  (4.33c)  we  first  introduce  the  following  change 
of  variaibles:  0 

f  =  X  -  Vt, 

y  =  y,  (5. la)  H 

ff  =  t. 

Derivatives  with  respect  to  x,  y,  t  can  be  expressed  as  ^ 

derivatives  in  the  new  variables  y,  a  as  follows: 


d  d 

■3^  " 

3  _  3 

■ly  “*  'y*  « 

dy  ^ 


d 

TE  " 


(5.  lb) 


Substituting  (5.1a,b)  into  (4.33c)  yields 


a^A, 


fJ!:+Vo,y,a) , 


(5.2) 
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where 


f.{S:+V(T,y,c)  =  c*  [  ^  ^ 

dr  dy^ 


(5.3) 


We  shall  solve  (5.2)  by  beginning  with  the  solution  to  the 
following  three-dimensional  wave  equation: 


_  2  r  d^U  d^U  d^u  1 

dt^  dx^  dy^  3z* 

(5.4a) 

u(x,y,  z,  0)  =<Kx,y,  z). 

(5.4b) 

• 

u,ix,y,2.Q)  ='^{x,y,  z). 

(5.4c) 

Kirchoff's  formula  (see  Strauss  1992), 

the  solution  to 

• 

c)  is 

iiix„,y„,2„,tj  =  ij—  rfi'(x,y,z)dS 

47rc^t„  ^ 

(5.5) 

• 
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where 


D  =  [  (x,y,z)  ;  (X-x,)*  +  (y-y^,)*  +  =  C^tl]  .  (5.6) 


Note  that  the  value  of  u(x„.y<,,  z,,,  t,)  depends  only  on  the  values  • 

of  ^{x,y,z)  euid  ^(x,y,z}  on  the  spherical  surface  D. 

By  the  method  of  descent,  we  now  obtain  the  solution  to 
the  following  two-dimensional  wave  equation:  • 


_  I ,  ^  d^u  , 

dt*  3x^  dy^  ' 

(5.7a) 

u(x,y,0)  =  *(x,y) , 

(5.7h) 

u,(x,y,0)  ='i'(x,y). 

(5.7c) 

We  regard  the  solution  u(x,y,t)  of  (5.7a-c)  as  a  solution  to 
the  three-dimensional  problem  (5.4a-c)  with  initial  data 
independent  of  z.  By  Kirchoff's  formula  (5.7),  we  have 


u(0,0,t)  =  — L_  r  r’J'{x,y)dS 
4irc*t  JpJ 


(5.8) 
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where 


A  =  nx.y.  2)  :  X*  +  y*  +  ar*  =  c*t*]  .  (5.9) 

This  is  the  solution  to  {5.7a-c)  but  it  can  be  8in5>lified.  To 
simplify  (5.8),  we  must  esqpress  the  surface  integral  as  an 
ordinary  double  integral  -  It  is  twice  the  integral  over  the 
top  heniisphere  z  »  (c^t^-x^-y^) where  (x,y)  ranges  over  the 
region  D,  -  {  (x,y)  :  x*  +  y^  s  c*t*  }.  By  definition 


dS  =  |J?,  X  J^|dx£fy, 


(5.10a) 


where 


and 


«  =  (  X,  y ,  /  { -y* )  ) , 


=  - - —  . 

>/(c^e-x2-y*) 


{5.10b) 


(5.10c) 


Substituting  (5.10a-c)  into  (5.8)  yields 


u(0,O,t)  .  .Jl-  f  f_ 

V  / 


^ix,y) 


cr*t*-x*-y* 


-.dxtfy 


Tt 


■X—  { 

r  . dxdv 

4, 

1  ^c^t^-x*-y* 

(5.11) 


This  is  the  solution  at  the  point  (0,0,t).  Therefore,  the 
solution  at  a  general  point  is 


1  r  f _ ^(3c,y) 


dxdy 


*(x,y) 


(5.12) 


y  ^cHl-  {x-xy-  {y-yy 


dxdy 


where 


D,  =  [  (x,y)  :  (x-xj^  +  (y-y„)^  *  (5.13)  ^ 

Note  that  (5.12)  shows  that  the  value  u(Xo,yo,to)  depends  on 
the  values  of  4>(x,y)  and  '♦(x,y)  inside  the  disc  defined  in  • 

(5.13).  This  illustrates  the  well  )aiown  fact  that  Huygens' 
principle  fails  in  two  dimensions.  We  shall  use  (5.12)  in  the 
solution  to  the  nonhomogeneous  two-dimensional  wave  equation  • 

discussed  in  the  next  section. 
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B.  mm»MXL*8  FBZHCXPLE 

We  now  consider  the  following  nonhcanogeneous  two- 
dimensional  wave  equation 


Hr  -  r  +  Hr  i 

3t*  dx*  dy^ 


gix.y,  t) , 


(5, 14a; 


with  homogeneous  initial  conditions 


u{x,y,0)  =  0, 

(5,14b) 

u,(x,y,0)  =  0, 

(5.14C) 

By  Duhamel's  principle,  the  solution  to  (5.l4a-c)  is  given  by 


u(x,y,  t)  =  JJ  £7(x,y,  t,6)  ds  (5.15) 


where  for  each  s  *  0,  U(x,y,t,s)  is  the  solution  of 


( 


3x2 


i!“)  .0 

dy^ 


for  ta  s. 


(5.16a) 


with  initial  data  prescribed  on  the  plane  t  -  s  such  that 
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U(x,y,a,3)  =  0, 


(5.l6jb) 


n,{x,y.3,s)  =  g(x,y,8)  .  (5,16c) 

Since  the  wave  operator  is  invariant  under  translations  in 
time  t  (see  Strauss  1992} ,  we  can  shift  the  time  t  by  a 
constant  and  still  have  a  solution.  We  let  V(x,y,t,s)  » 
U(x,y, t+s, s) ,  so  that  V(x,y,t,s)  is  a  solution  to 


V„-C^{V^*V„)  =0  for  t  a  0.  (5.17a) 


with  initial  data  prescribed  on  the  plane  t  =  0  such  that 


V(x.y,0,s)  *  n{x,y,s.8)  =  0.  (5.17b)  ^  ^ 

V,(x,y,0,s)  =  V,{x,y,s,3)  =  g(x,y,3)  .  (5.17c) 

Using  Kirchoff's  formula  (5.12)  to  obtain  a  solution  to 
{5.l7a-c)  then  yields 

V(x,y,t,8)  =  -i-  r  f-  ■  .  dxdv  (5.18) 

T>,  Vc*t^-  (x-x)*-  (y-y)  * 
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where  D,  =  {  {x,y)  ;  (x-x)^  +  (y-y)^  *  cV  ).  Since  our  purpose 
here  is  to  solve  the  nonhc»itogeneous  two-dimensional  %rave 
equation  in  {5.14a-c)  by  Duhamel's  principle  (5. 15),  we  can 
again  shift  the  time  t  by  a  constant  so  that  U(x,y,t,s)  « 
V(x,y, t-s,s) .  This  last  expression  for  U(x,y,t,a)  together 
with  (5.18)  then  inplies 


U{x,y,  t,s) 


1  f  f _ g(x.y.s) _ 

2irc  ^ cr2  ( t-s)  (x-x)  (y-y)  ^ 


dxdy 


(5.19) 


where 


A  t  (x,y)  :  (x-x)^  +  iy-y)^  s  c^it-s)^  ]  .  (5.20) 


Now  using  Duhcunel's  principle  (5.15),  we  obtain 


•  • 


u{x,y,t)  = 


2irc 


UJ 


g{x,y,  8) 


<}  ( t-s)  (x-x)  iy-y)  ^ 


dxdy 


ds  (5.21) 


the  solution  to  (5.14a-c)  where  D,  is  defined  in  (5.20). 

In  deriving  the  solution  to  the  nonhomogeneous  two- 
dimensional  wave  equation (5. 14a-c) ,  we  have  solved  (5.2).  We 
change  the  names  of  the  dependent  and  independent  variedsles  in 
(5.21)  to  obtain  the  solution  of  (5.2): 
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da  (5.22) 


A,{t.y,o) 


^  f 

V  r  .  .  ..  csd?] 

2irc4 

where  D,  »  {  (x,y)  :  (x-f)*  +  (y-y)*  *  c^(ff-s)^  }-  Letting  - 
X  -  Vt  where  a  «  t,  (5.22)  bec<xnes 


A,{x.y,t)  = 

41/ 


2irc. 


fJx+V8,y.  8) 


( t-s)  (x-  (x-vt) )  ^-  (y  y)  ^ 


dxdy 


ds 


[5.23) 


where  D,  =  {  {x,y)  :  (x-{x-Vt))^  +  (y-y)^  s  c,^{t-s)^  }. 

To  derive  the  solution  for  the  nonhomogeneous  two- 
dimensional  wave  equation,  we  began  by  solving  the  homogeneous 
three-dimensional  wave  equation  using  Kirchoff's  formula 
(5.5)  .  Next,  we  used  the  method  of  descent  to  solve  the 
homogeneous  two-dimensional  wave  equation  by  reducing 
Kirchoff's  surface  integral  formula  to  an  ordinary  double 
integral  (5.12) .  We  developed  the  solution  to  the  homogeneous 
two-dimensional  wave  equation  in  order  to  solve  the 
nonhomogeneous  two-dimensional  wave  equation  by  applying 
Duhamel ' s  principle  (5.21).  This  gave  the  solution  of  the 
wave  equation  in  (5.2)  in  the  form  given  above  in  (5.23)  .  In 
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Chapter  VI  we  will  examine  the  solution  to  {5.23}  for  long 


VZ.  LONG  TIKE  SOLGTZC»7 


A.  C,^  7t 

Before  looJcing  at  the  behavior  of  (5.23)  as  t  -*  »  (where 
the  topography  f, (x+Vs,y,s)  has  a  general  form),  we  shall 
first  introduce  a  specific  moving  point  source  for  the 
topography,  so  that 

f,  ix,  y,B)  =6  (x-x,)  6  (y-y„)  g  ( s)  . 

Consequently 

f,(x+ys,y,  s)  =  6(x+V’s-xJ  5(y-y„)  g{s) , 

where  g  is  a  non-negative  function  satisfying 

g(0)  =  0  and  g(a)  =  1  for  s  a  s„.  (6.1c) 

If  we  can  solve  (5.23)  by  substituting  (6.1b)  in  for 
f, (x+Vs, y,  s)  ,  then  we  can  solve  (5.23)  for  any  form  of 
topography  by  applying  the  principle  of  superposition. 

To  solve  (5.23),  we  begin  by  substituting  (6.1b)  into 
(5.23)  to  obtain 


(6.1a) 


( 6 .  lb) 
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A^(x,y,t)  = 


r  r  6  (jc+Vs-xJ  a  (y-y^)  g(s) 
y  Jcj ( t-s)  [X-  {X- Vt)  ]  (y-y)  ^ 


dxdy  ds. 


(6.2) 


We  let  X  -  (x+Vs-x„)  and  Y  •  (y-Yo)  and  substitute  into  (6.2) 
to  obtain 


A,ix,y,  t) 


r  d{X)  6{Y)  g-(g)  dXdY  ^ 

c]  ( t-s)  [X+V(  t-s)  -  (x-x^)  ]  ^-  [  Y-  {y-yj  ]  ^ 


(6.3) 


where 


D,  *[(X.Y)  :  [X+V(C-S) -(X-XJ3’  +  *  c/(t-s)*].  (6.4) 


In  order  to  solve  (6,3)  we  must  determine  those  s  for  which 
(6.4)  is  satisfied.  Letting  t  »  (t-s)  for  (t-s)  2  0,  (6.4) 
sinplifies  to 


(c/-v*)t*  +  2V{x-x„)7  -  [(x-x„)*  +  (y-y„)*]  »  0.  (6.5) 
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Note  that  X  and  Y  are  not  included  in  (6.5)  since  the 
integration  in  (6.3)  with  respect  to  X  and  Y  makes  both  X  and 
Y  become  zero.  Solving  (6.5)  for  t  yields 


Tj.  = 


-V(x-xJ  ± 


+  (c;-v®)  (y-yj‘ 


c.-V^ 


(6.6) 


The  values  of  and  r.  will  depend  on  whether  c,*  >  or  c,^  < 
V^.  We  shall  exounine  both  cases  in  the  following  sections. 

1.  C.^  > 

When  the  long -wave  phase  speed  is  greater  than  the 
velocity  of  the  submarine,  we  see  from  (6.6)  that  t+  >  0  and 
T.  <  0.  Factoring  (6.5)  yields 

icJ-V^)  (T-T.)  (T-Tj  *  0  (6.7) 


where  t+  and  t.  are  defined  in  (6.6)  .  Since  the  left  hand  side 
of  (6.7)  must  be  greater  than  or  equal  to  zero,  t  must  be 
greater  than  or  equal  to  (r  %  r+) .  Hence,  s  must  be  less 
than  or  equal  to  (t-T+).  Therefore,  the  upper  limit  of  s  in 
(6.3)  becomes  (t-r+).  Now,  integrating  (6.3)  with  respect  to 
X  and  Y  yields 


51 


(6.8) 


A.{x,y,t)  =  — \ — r\ 

2xc,(c,  V*)  ®  /(( t-r*)  -s)  ({ t-T  J  -s) 


where  t+  and  t.  are  defined  in  (6.6).  Using  (6.1c),  we  break 
(6.8)  into  the  following  parts 


A,(x,y,  t)  =  J,  +  Jj, 


(6.9a) 


where 


j  ^  _ 1  _ g(a)  ds 

2irc,(c,-V^)''^  “  ^((t-Tj  -s)((  t-1 


.)-fl) 


(6.9b) 


and 


I,  * 


i: _ f'-'* 


ds 


2irc,(c,-l^)‘^*>.  /((t-Tj-s)((t-T.)-s) 


(6.9c) 


We  first  look  at  the  behavior  of  I,  as  t  -*  ».  We  see 
from  (6.9b)  that  for  long  time,  I,  is  0(l/t)  .  Therefore,  we 
can  neglect  this  portion  of  A,(x,y,t).  Now  in  evaluating  • 

for  long  time,  we  let  r  <■  (t-s)  so  that  Ij  simplifies  to 
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I  =  3.  n 


« dr 


^(r-T4(T-r.) 


(6.10) 


We  cc»nplete  the  square  of  the  denominator  of  (6.10)  to  obtain 


h  =  - 1 - f‘**'-7= 


C*T 


(6.11) 


■a)’-Jb’ 


Here  a  -  (t++t.)/2,  and  -  (t  +  -t.)V4.  Letting  bu  -  (r-a)  and 
substituting  into  (6.11)  then  yields 


^2 


_ 1 _ 

2irc,(c,*-V2)‘'* 


du 

Vu^-l 


(6.12) 


Integrating  Ij,  we  find 


1 

,„i  . 

•  ' 

(t-s,-a)*  ^ 

i;2 

1 

2irc.{c,^-V®)‘'^ 

1  jy 

b* 

1 

+ 


■1 


• 

i;z 

Inl 

{r.-a)*  , 

1 

'  b 

b* 

) 

(6.13) 


:  • 


For  long  time  t,  Ij  is  0(ln  t) . 

In  the  case  where  c,*  >  V^,  we  see  that  the  long  time 
solution  for  A,(x,y,t)  is  0(ln  t)  .  Next  we  examine  the  cose 
for  which  c,^  <  V*. 

2.  C.*  <  V* 

When  the  long-wave  phase  speed  is  less  than  the 
velocity  of  the  submarine,  we  see  from  (6.5)  that  x  must  be 
greater  that  x^  since  t  is  greater  than  or  equal  to  zero.  We 
must  reexamine  (6.6)  to  see  if  a  positive  real  value  for  r  is 
possible.  Rearranging  (6.6)  for  c,^  <  V*  yields 


v(x-x„)  ±  ^c/{x-x^)^  -  (V^-cJ)  (y-y„)’ 
V^-cJ 


(6.14) 
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Prom  (6.14)  we  see  that  a  real  solution  for  t  exists  only  if 
the  following  inequality  holds 


Since  T+  and  r.  are  both  positive,  we  see  from  (6.7)  that  t 
must  be  between  and  t.,  ie.,  (t.st*t+)  .  Therefore,  s  must 
be  less  than  or  equal  to  (t-r.)  • 

To  solve  A,(x,y,t)  for  c,^  <  as  t  we  start  with 
(6.3).  Integrating  (6.3)  with  respect  to  X  and  Y  yields 
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A^ix,y,t)  = 


I 


gia)  da 


2irc,{V®  Ct)  -6  fir)({t-r.)  -a) 


<6.17) 


where  t+  suid  t.  are  defined  in  (6.14).  Using  (6.1c),  we  break 
(6.17)  into  the  following  parts 


A,  (x,  y,  t)  =  Jj  +  J4  , 


(6.18a) 


where 


= 


1  f 


_ . _ gia)  da _ 

2ircJXr^~c,)^'^i  -e)((t-Tj  -a) 


(6.18b) 


and 


= 


«  T 

_ 1 _  r _ da 

2irc,{c,-V^)^'^  I  ^((t-T,)  -s){(t-Tj  -s) 


(6.18c) 


Looking  at  the  behavior  of  I3  as  t  ®,  we  see  from 
(6.10b)  that  I3  is  0(l/t).  Therefore,  we  can  neglect  this 
portion  of  A,{x,y,t).  The  evaluation  of  I4  is  similar  to  the 
analysis  as  developed  for  Ij: 
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1/2 

1 

(t-s^-a) 

Ini  "  ^ 

1 

tu  1  ^ 

. 

1 

+ 


-1 

2T(c,{v^-c^,y'^ 


-  1  (T.-a) 

1  n  I 

(T.-a)^  , 

Ill 

1 

'  b 

I 

(6.19) 


Here  a  »  <t^.+t.)/2  and  -  (t+-t.)^/a.  For  long  time,  I,  is 
0(ln  t)  . 

In  the  case  where  V’  >  c^,  we  see  that  the  long  time 
solution  for  A,{x,y,t)  is  0(ln  t)  .  Next,  we  examine  the 
resonant  case  where  c,^  *  V^. 

B.  C.*  -  V* 

Previously,  we  developed  the  long  time  solution  for  the 
cases  when  c,*  >  and  c,^  <  with  the  interesting  result  that 
in  both  cases  the  solution  had  growth  of  0(ln  t) .  The  reason 
for  this  growth  is  the  failure  of  Huygens'  principle  in  two 
dimensions.  Though  ln(t)  tends  to  infinity  for  long  time,  it 
is  a  function  of  extremely  slow  growth,  and  for  the  times  we 
are  interested  in,  it  can  be  considered  as  a  bounded  function. 

In  the  resonant  case,  we  are  concerned  with  the  behavior 
of  A,(x,y,t)  as  t  -*  00  for  the  case  when  the  speed  of  the 
s\abmarine  is  equal  to  the  long-wave  phase  ^peed  (c,*  =  V^)  . 


57 


To  determine  if  a  solution  to  (6.3)  is  possible  for  c,*  = 
V*,  we  naist  reexamine  (6.5).  From  (6.5)  we  see  that  a 
solution  exists  only  if  s  satisfies  the  following  inequality 


s  *  t  - 


(x-x^)2  +  (y-yj 
2V(x-jr) 


for 


V(x-x-)  >  0. 


(6.20) 


Integrating  (6.3)  with  respect  to  X  and  Y  yields 


A,(x,y,  t)  = 

1  f'*  g(s)  ds  (6.21) 

2rcJo  I 

pvix-xa)  it-s)  -  i(x-xy  ^  {y-y,)M 


Here  A  »  [t  -  ( {x-x„)’+ (y-y,)’)  /  {2V(x-Xo) )  ]  .  Using  (6.1c),  we 
break  (6.21)  into  the  following  parts 

Ajx.y.  t)  =  Jj  +  Ij  ,  (6.22a) 


where 


X,  = 

1  r _ g(g)  dg _ 

^  \l2V(x-x^)  {t-B)  -  [{x-x„)2+  (y-yj^l 

and 


(6.22Jb) 


^6  = 

_ ds _ <"*220 

*.  ^2V{x-x^)  It-s)  -  [(x-x„)^  +  <y-y;,)^  l 

Here  A  =>  [t  -  ( (x-Xo)^+ (y-Yo^^)  /  <2V(x-x„) ) }  . 

Looking  at  the  behavior  of  I5  as  t  -*  »,  we  see  from 
{6,22b)  that  for  long  time,  Ij  is  Od/t*'*)  .  Therefore  we  can 
neglect  this  portion  of  A,{x,y,t)  .  Integrating  Ij  directly,  we 
obtain 


^6  = 

I - ; - (6.23) 

1  ]]2V[x-xj  [t-sj  -  [{x-xj^  +  (y-y„)^] 

2  ire,  V(x-x^) 


From  (6.23)  we  see  that  Ij  is  0(t''*)  as  t  -»  00.  Therefore,  we 
see  that  the  solution  for  A,(x,y,t)  at  resonance  grows  like 
©{t''^)  for  long  time. 


m 


m 
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To  derive  the  long  time  solution  (t  -•  «)  for  A,{x,y,t),  we 
began  by  introducing  a  moving  point  source  for  the  topography, 
f,(x+Vs,y,s) ,  with  the  knowledge  that  if  a  solution  was 
attainable  by  using  {6.1c),  we  could  solve  (5.23)  for  any  form 
of  the  topography  by  the  principle  of  superposition.  We  first 
looked  at  the  case  where  c,^  >  V*  and  determined  that  the  long 
time  solution  for  A.{x,y,t}  was  0(ln  t)  .  Next,  we  looked  at 
the  case  where  c,^  <  euid  saw  that  a  solution  existed  only 
within  a  certain  region  defined  by  (6.16a,b) .  The  long  time 
solution  for  A,(x,y,t)  in  this  case  was  also  0(ln  t)  .  We 
attribute  the  0{ln  t)  growth  to  the  failure  of  Huygens' 
principle  in  two  dimensions.  For  the  resonant  case  where  c,^ 
V^,  we  found  the  solution  for  A,{x,y,t)  to  be  OCt*'^)  for  long 
time.  From  this  analysis,  we  see  in  this  problem  that  if  a 
submarine  patrols  at  or  near  a  long-wave  phase  speed  a 
significant  resonant  could  occur  for  long  times. 


VIX.  CONCLUSIONS 


This  thesis  has  presented  a  theoretical  study  of  the  flow 
of  a  stratified  fluid  over  variable  topography  as  a  loodel  of 
the  resonant  interaction  of  a  sulmarine's  waUce  with  a 
stratified  fluid.  We  considered  the  case  in  which  a  subotarine 
patrols  at  a  speed  near  the  long-wave  phase  speed  of  an 
internal  wave.  Internal  waves  are  generated  by  the  flow 
around  the  submarine’s  wake. 

The  long  time  solution  obtained  in  our  model  for  the  cases 
when  c,^  >  and  had  gro»#th  of  0(ln  t)  .  (see  (6.13)  and 
(6.19)).  We  attribute  the  0(ln  t)  growth  to  the  failure  of 
Huygens*  principle  in  two  dimensions.  For  the  resonant  case 
where  c,^  =  V^,  we  found  the  solution  to  be  0(t''^)  for  long 
time.  (see  (6.23)).  It  is  important  to  note  that  for  the 
time  scales  we  are  concerned  with,  both  long  time  solutions 
for  the  nonresonant  cases  are  bounded  and  have  extremely  slow 
growth.  Therefore  it  is  the  resonant  case  which  is  important. 

The  growth  in  time  occurring  in  the  resonant  case 
indicates  that  perhaps  nonlinear  effects,  balanced  by 
dispersion,  could  cause  significcmt  upstream  disturbances. 
Frcm  these  disturbances,  the  possibility  exists  of  detecting 
a  sulxnarine  by  the  internal  waves  its  wake  generates. 
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Further  studies  in  this  problen  should  begin  by  rescaling 
the  asyaqjtotic  esqpansions  to  balance  the  leading  order 
quadratic  nonlinear  terms  with  the  topographic  forcing.  This 
will  achieve  a  balance  between  nonlinearity  and  dispersion. 
We  expect  the  solution  to  this  problem  to  show  that 
significant  upstream  disturbcuices  could  occur  when  the 
su]ixnarine  patrols  at  a  speed  at  or  near  a  long-wave  phase 
speed. 

If  internal  %#aves  were  generated  only  by  underwater  moving 
objects  then  the  problem  of  detecting  the  submarine  by  the 
internal  waves  its  wake  generates  would  be  made  simpler,  but 
this  is  not  the  case.  Internal  %»aves  are  a  natural  occurring 
phencanenon  in  the  oceems.  This  complicates  the  problem 
because  in  order  to  detect  the  siatanarine  by  the  internal  waves 
that  it  generates,  the  ambient  internal  wave  patterns  at  the 
location  of  the  submarine  must  be  known.  The  problem  of 
discerning  ambient  internal  waves  frcan  internal  waves 
generated  by  the  submarine's  wake  is  similar  to  the  problem  of 
detecting  a  submarine  passively  (such  as  with  a  tail) .  The 
sonar  operator  filters  out  the  ambient  noises  or  background 
noises  of  the  ocean  from  the  incoming  signal  to  detect  the 
acoustic  signature  of  the  submarine.  From  further 
understemding  of  the  ocean's  naturally  occurring  internal 
waves,  we  may  one  day  be  cdsle  to  discern  internal  waves 
generated  by  a  s\ibmarine  frcxn  those  of  the  ocean. 
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